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Abstract 

As an alternative to dark energy that explains the observed acceleration of the universe, it has 
been suggested that we may be at the center of an inhomogeneous isotropic universe described by 
a Lemaitre-Tolman-Bondi (LTB) solution of Einstein's field equations. To test this possibility, it is 
necessary to solve the null geodesies. In this paper, we derive a fully analytical set of differential 
equations for the radial geodesies as functions of the redshift in LTB models. Then we use them to 
show that a positive averaged acceleration ao obtained in LTB models through spatial averaging can 
be incompatible with cosmological observations. We provide examples of LTB models with positive 
a£) which fail to reproduce the observed luminosity distance Dl{z). Since the cosmic acceleration 
^FLRW jg obtained from fitting the observed luminosity distance to a FLRW model we conclude that 
in general a positive a/) in LTB models does not imply a positive a^^^^ . 
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I. INTRODUCTION 



High redshift luminosity distance measurements 2l|, |28!, |29|, l30|, l3l|, l32] and the WMAP 



measurement 
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35| of cosmic microwave background (CMB) interpreted in the context of 



standard FLRW cosmological models have strongly disfavored a matter dominated universe, and 
strongly supported a dominant dark energy component, giving rise to a positive cosmological 
acceleration, which we will denote by a^^^^ (not to be confused with the scale factor a). As 
an alternative to dark energy, it has been proposed that we may be at the center of 

an inhomogeneous isotropic universe described by a Lemaitre-Tolman-Bondi (LTB) solution 
of Einstein's field equations, where spatial averaging over one expanding and one contracting 
region is producing a positive averaged acceleration Another more general approach to 
map luminosity distance as a function of redshift Dl{z) to LTB models has been recently 



proposed [23 



2J], showing that an inversion method can be applied successfully to reproduce 



the observed Dl{z). 

The main point is that the luminosity distance is in general sensitive to the geometry of 
the space through which photons are propagating along null geodesies, and therefore arranging 
appropriately the geometry of a given cosmological model it is possible to reproduce a given 
Dl{z). For FLRW models this corresponds to the determination of and and for LTB 
models it allows to determine the functions E{r), M{r),tb{r). 

The averaged acceleration a^) on the other side is not directly related to a^^^^ , since the 
latter is obtained by integrating the position dependent cosmological redshift along the null 



geodesies, while a/j is the result of spatia. 
of the underlying space. It was shown 



averaging, and has no relation to the causal structure 
25j that this crucial difference can make a £, unobservable 
to a central observer Oc when the scale of the spatial averaging is greater than its event horizon. 

In this paper we will further investigate the relation between LTB models with positive 
averaged acceleration and cosmological observations, showing that in general they can be in- 
compatible. 



II. LEMAITRE-TOLMAN-BONDI (LTB) SOLUTION 

Lemaitre-Tolman-Bondi solution can be written as 



-dt 



2 , CM!^^2^^2 

1 + 2E 
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where R is a function of the time coordinate t and the radial coordinate r, R = R{t,r), E is 
an arbitrary function of r, E = E{r) and R,r = OR/ dr. 
Einstein's equations give 

' 2E{r) ^ 2M(r) ^ 



R. 



p{t,r) 



i?2 

R^Rjr 
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(3) 



with M = M{r) being an arbitrary function of r and the dot denoting the partial derivative 
with respect to t, R = dR{t, r)/dt. The solution of Eq. ([2]) can be expressed parametrically in 
terms of a time variable i] = dt' /R{t', r) as 

M(r) 



R{r],r) 



-2E{r) 

M{r) 
-2E(r) 



1 - cos ( ^-2E{r)r] 
1 



T] 



-2E(r) 



sin J —2E{r)ri 



(4) 
(5) 



where R has been introduced to make clear the distinction between the two functions R{t, r) 
and R{ri, r) which are trivially related by 



R{t,r) = R{r]{t,r),r) 



(6) 



and tfc(r) is another arbitrary function of r, called the bang function, which corresponds to the 
fact that big-bang/ crunches can happen at different times. This inhomogeneity of the location 
of the singularities is one of the origins of the possible causal separation 25| between the central 
observer and the spatially averaged region for models with positive ao- 
We introduce the variables 

2E(r) , . 6M(r) 



Po[r) 



(7) 



SO that Eq. ([1]) and the Einstein equations ([2]) and ([3]) are written in a form similar to those 
for FLRW models, 
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The solution of Eqs. (jl]) and 

a(?7,r) = 



]) can now be written as 
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where f] = r]r = dt' / a{t' , r) . 

In the rest of paper we will use this last set of equations and drop the tilde to make the 
notation simpler. Furthermore, without loss of generality, we may set the function poir) to be 
a constant, po{r) = po = constant. 



III. GEODESIC EQUATIONS 

The luminosity distance for a central observer in a LTB space as a function of the redshift 
is expressed as 

Dl{z) = (1 + zfR (t(z), r{z)) = (1 + zfr{z)a {v{z), r{z)) , (13) 

where (t{z),r{z)^ or (j_^ri{z),r{z)^ is the solution of the radial geodesic equation as a function 
of the redshift. The past-directed radial null geodesic is given by 

dTjr) -Rr{t,r) 

— — = / T r ,r ; / t,r = ^ 14 
dr + 2E{r) 

Using this equation and the definition of redshift, we can derive the geodesic equations for the 
coordinates (?], r), 

drj{z) dt{ri,r)/dr — f 



dz (1 + z)af 

dr{z) 1 



(15) 
(16) 



dz {l + z)f 

where / = df{t,r)/dt. 

The derivation of the analytical solution is based on the use of the conformal time variable 
rj, which by construction satisfies the relation, 

= a-' . (17) 
dt ^ ' 

This means 

t{r],r) = tb{r) + / a{r] ,r)dT] , (18) 
Jo 

dt = a{ri, r)dri + ( ^"^^ '^^ c?V + tj,(r) ) dr , (19) 







dr 



In order to use the analytical solution we need to express / and / in terms of rj. This can 
always be done by using 



t{0,r)=h{r), a(0,r)=0 
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^ drjdr drf dr drj drj \dr J J 

In this way the coefficients of equations f|T5|) and f|T6|) are fully analytical, which is a significant 
improvement over previous approaches, and could be used to give analytical solutions in the 
low redshift regime. 



IV. AVERAGED ACCELERATION 

Following the standard averaging procedure [3, 126I, |27| we define the volume for a spherical 
domain, < r < re, as 

Vo = 4vr/ dr, (26) 

^0 Jl + 2E(r) 



and the length associated with the domain as 

Ld = (27) 

Then the deceleration parameter and the averaged acceleration a£, (not to be confused with 
the scale factor a) are defined as 

qD = -LdLd/LI, (28) 
an = Ld/Ld. (29) 

As models that give a positive averaged acceleration, we consider those studied in p^]. They 
are characterized by the functions k{r) and tb{r) given by 

^ ^ 1 + (r/r^)"* ^ ^ 

^^'^ " 1 + (r/r,)- + ^ ■ ^^^> 



Note that they fix the length scale by setting k{0) = 1. After exploring the 9 parameters space 
they give three examples of LTB solutions with positive as shown in Table I. 



TABLE I: Three examples of the domain acceleration. 
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Defining tq as the time in the first column of Table I, i.e. the time at which q{tq) = qd, with 
qd being the value in the last column of the same table, we solved the null geodesies equations 
(ITSll and (JT6ll imposing the following initial conditions. 

r]iz = 0) = r/, = r/(t„0), (32) 
r{z = 0) = 0, (33) 

where ri(tq,0) is obtained by solving numerically for rjg, Eq. (IT^ . 

t{Vg,0)=tq, (34) 

For the second and third model of Table I which we will study we obtain, respectively, 

r]g = 0.0330199, (35) 
% = 3.30189 X 10"^ (36) 

Eq. (l32l) is the natural way to map these models into luminosity distance observations for 
a central observer which should receive the light rays at the time tq at which the averaged 
acceleration is positive. 



V. COMPARING LTB TO OBSERVATIONS 



An important model independent quantity is the expansion history, H{z), whose value can 
be reconstructed from observations of the luminosity distance, D^, via a single differentiation 
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FIG. 2: k{r{z)) and p{rj{z),r{z)) are plotted for the second model in Table I. 



Note that this assumes the universe is not only homogeneous and isotropic but also spatially 
fiat. The equation of state, w{z), of DE is more cumbersome to reconstruct since it involves 
second derivatives of Dl{z) and is therefore a noisier quantity than H{z). An additional source 
of uncertainty relating to w{z) is caused by the fact that the value of the matter density, Qom 
enters into the determination of w{z) explicitly, through the expression 

, , (2x/3) d InH / dx — 1 , , 

'"-(^) - 1 - (H^/nk^ .3 ■ (38) 
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FIG. 3: R{r]{z),r{z)) and R^^.{r]{z),r{z)) are plotted for the second model in Table I. 
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FIG. 4: H{z) and Di{z) are plotted for the second model in Table I. The dashed line corresponds to 
the concordance ACDM model. 

In our case we can first calculate Di{z) according to the LTB solutions geodesies equation we 
derived in the previous section, and then use it construct H{z) and w{z) according to Eqs. fl38|) 
and (1371). 

We compare H{z) and Dl{z) to the observed Hp and Df^ using the standard best fitted 
FLRW cosmological model, 

Hpizf = Hl(n^{i + zf + r]„ = o.25, (39) 

Dl{z) = {l + z) r -^dx, (40) 
It should be noted that the parameters given in Table I require to fix the units, and we do that 
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z 



FIG. 5: ■w{z) and q{z) are plotted for the second model in Table I. 




z 



FIG. 6: wde{z) and a{z) are plotted for the second model in Table I. 



by imposing 

H{0)=Ho. (41) 

In this way we eliminate any ambiguity in the choice of units in order to be able to compare 
the observed Hf{z) and H{z). 

For both the second and third models in Table I, H{z) and Dl{z) are not reproducing 
correctly observational data both from the quantitative and qualitative points of view. The 
effective w{z) is small and positive and the effective W£ie{z) is not —1 for both the second and 
third models. On the contrary, for the redshift range of z ^ 1, wde{z) > 1 for the second 
model, while 1 ^ wde{z) > for the third model. 

As it can be seen from the plots, k{z) is approximatively constant in the observationally 
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FIG. 7: tf){r) and k(r) are plotted for the third model in Table I. 
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FIG. 8: k{r{z)) and p{rj{z),r{z)) are plotted for the third model in Table I. 



interesting redshift range, but given the parameters of these models matter is dominating at 
the time tg, since as shown above rjq <^ 1, which explains why w is so small and q{z) ^ 1/2. 



z) com- 



These examples show how a positive does not imply a luminosity distance Dr\ 
patible with observations, and gives a reverse example of the results obtained in ll], where 
they obtained a LTB model which fits the observed luminosity distance, consistent with a pos- 
itive a^^^^ , but without positive averaged acceleration ao. Our results do not rule out LTB 
models as alternatives to dark energy since the inversion method 23|, |2J| allows us to obtain 



the observed luminosity distance without any averaging, and some concrete examples derived 



independently have already been proposed 
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FIG. 9: R{r]{z),r{z)) and R^r{ri{z),r{z)) are plotted for the third model in Table I. 
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FIG. 10: H(z) and Dl{z) are plotted for the third model in Table I. The dashed line corresponds to 
the concordance ACDM model. 



VI. DISCUSSION 

We have derived a set of differential equations for the radial null geodesies in LTB space- 
time without a cosmological constant and used them to compute the luminosity distance for 
models which have positive spatial acceleration. We have also extended our analysis to other 
observables such as H{z) and w{z), explicitly showing that these models are not compatible 
with observations. Our conclusions do not rule out LTB models as alternatives to dark energy, 
but provide further evidence that physical quantities obtained via spatial averaging are not 
relevant to explain observational data. 

In the future it will be interesting to find general analytical results in support of our con- 
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FIG. 11: w{z) and q{z) are plotted for the third model in Table I. 
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FIG. 12: wde{z) and a{z) are plotted for the third model in Table I. 



elusions and to use the equations we derived to provide a new method to solve the inversion 
problem of mapping the observed luminosity distance Dl{z) to LTB models. 
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